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Abstract 

We obtain gravity duals to supersymmetric gauge theories in two and three spacetime di- 
mensions with unquenched flavor. The supergravity solutions are generated by a set of color 
branes wrapping a compact cycle in a Calabi-Yau threefold, together with another set of 
flavor branes extended along the directions orthogonal to the cycle wrapped by the color 
branes. We construct supergravity backgrounds which include the backreaction induced by 
a smeared set of flavor branes, which act as delocalized dynamical sources of the different 
supergravity flelds. 
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1 Introduction 

The AdS/CFT correspondence [H [21 [3l H] is one of the major achievements of theoretical 
physics of the last few years. Indeed, this correspondence has provided analytic tools to 
explore the strong coupling regime of a large variety of gauge theories in the planar limit 
Ac — > oo. Originally, the AdS/CFT correspondence was formulated as a duahty between 
A/" = 4 super Yang-Mills (SYM) gauge theories in four-dimensions and type IIB supergravity 
in AdS^ X S^. Extending the duality beyond this highly symmetric theory is clearly a 
problem which deserves to be studied and, for this reason, there has been much effort devoted 
to constructing supergravity duals to theories with lower amounts of supersymmetry and 
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different field content in several spacetime dimensions. 

A general approach to construct gravity duals with lower amounts of supersymmetry 
consists of considering higher dimensional branes wrapping cycles inside a non-compact 
manifold of special holonomy [3]. At energies small compared with the size of the cycles 
these solutions provide gravity duals of SYM theories living on the unwrapped dimensions 
of the brane. A very useful strategy to find supergravity solutions corresponding to wrapped 
branes is the use of the appropriate gauged supergravities, which are lower-dimensional 
theories in which the brane behaves as a domain wall. In these theories one can identify the 
spin connection along the wrapped cycle with a gauge field, implementing in this way the so- 
called topological twist [6]. In this approach the ten-dimensional background is obtained from 
the lower dimensional one by using the corresponding uplifting formulae. This program has 
been successfully carried out to find duals to gauge theories in several spacetime dimensions 
with different amounts of supersymmetry [7]-|17]. 

Another important generalization of the AdS/CFT correspondence has been the addition 
of matter degrees of freedom transforming in the fundamental representation of the gauge 
group, i.e. quarks. This new matter sector can be added by including flavor branes, which 
extend along all gauge theory directions and, in order to make its worldvolume symmetry a 
global symmetry from the gauge theory point of view, they should be also extended along 
some other non-compact directions [18]. When the number Nf of flavor branes is small 
compared with the number Nc of colors, one can treat the flavor branes as probes which do 
not modify the background created by the color branes. This defines the so-called quenched 
approximation which, on the field theory side, corresponds to considering the quarks as 
external fields that do not run in loops and, thus, to neglect the quantum effects produced 
by the fundamentals. By studying the worldvolume physics in this probe approximation one 
can address many interesting problems such as, for example, the meson spectrum [19] (see 
[20] for a review). 

When the number of flavors is of the order of the number of colors {Nf ~ Nc) one 
must face the problem of studying the deformation produced by the flavor branes on the 
geometry. Computing this backreaction is, in general, a very complicated task. Indeed, 
in these unquenched setups the flavor branes should be regarded as dynamical sources of 
the different supergravity fields and, unless one adopts some simplifying assumptions, the 
problem becomes, in many cases, intractable. For this reason, in this paper we will consider 
the case in which the flavor branes are homogeneously distributed by forming a smeared set. 
This smearing technique to add unquenched flavor was first introduced in [21] in a non-critical 
string framework and in [22J in a well-controlled ten-dimensional context. Subsequently, this 
approach has been successfully applied in several brane setups [23]- [33] (see [M] for a review). 

In this paper we study the addition of unquenched flavor to backgrounds constructed by 
wrapping D3-, D4- and D5-branes along two- and three-cycles of a Calabi-Yau threefold. 
The unflavored backgrounds, which are dual to gauge theories in two and three dimensions. 
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have been constructed in refs. [5l [HI [15] from the appropriate gauged supergravity. Here 
we reformulate these models in terms of a new set of variables which greatly simplify their 
interpretation and allows the introduction of the corresponding flavor deformation in a neat 
way. After this reformulation the background is determined by a set of functions depending 
on two variables, which represent the radial variables inside the Calabi-Yau cone and in the 
space transverse to it. These functions satisfy a set of first-order differential equations, which 
can be recast compactly in terms of a generalized calibration form. 

The flavor branes for the setups considered here are extended across the directions normal 
to the cycle that the color branes wrap. In order to obtain the backreaction of these flavor 
branes on the background the key point is realizing that they act as sources for the RR field 
and, as a consequence, they induce a violation of the corresponding Bianchi identity. For a set 
of delocalized flavor branes the violation of the Bianchi identity compatible with the preserved 
supersymmetry can be encoded in a simple modification (parametrized by two functions) of 
the ansatz of the RR form, while the ten-dimensional metric has the same form as in the 
unflavored case. By requiring certain regularity conditions, the functions which parametrize 
the flavor deformation can be greatly constrained and, in fact, they can be exactly matched 
with those obtained by a careful microscopic counting of the RR charge distribution produced 
by a family of supersymmetric embeddings. In this paper we will determine these embeddings 
and we shall perform the matching with the macroscopic analysis for the three cases analyzed. 
Moreover, the backreacted backgrounds will be found by a numerical integration of a PDE 
with sources. 

The rest of this paper is organized as follows. In section [2] we will consider our first 
background, namely the one generated by D3-branes wrapping a two-cycle in a Calabi-Yau 
threefold and preserving four supercharges. This background is dual to a two-dimensional 
gauge theory with = (2, 2) supersymmetry. We first reformulate the results of [5] for 
the unflavored case by introducing a second-order master differential equation whose solu- 
tions determine the unflavored background in implicit form. Then, we analyze the flavor 
deformation induced by a smeared set of D3-branes extended along two of the directions 
normal to the cycle within the Calabi-Yau cone. This deformation is analyzed first from 
the macroscopic point of view and, then, these results are reproduced from a microscopic 
analysis of the embeddings. Section [2] ends with the presentation of the numerical results 
for the flavored model. 

Section [3] is devoted to the study of a gravity dual to A/" = 2 gauge theories in three- 
dimensions generated by wrapping D4-branes on a two-cycle of a Calabi-Yau cone of com- 
plex dimension three. The analysis in this section runs completely in parallel with the one 
performed in section [2j In section [H we consider the backgrounds with four supersymmetries 
which are generated by D5-branes wrapping a three-cycle. The corresponding unflavored 
solution was found in refs. [m [15], while the macroscopic approach to adding flavor was 
first introduced in ref. [30]. In section [D we identify the deformation which corresponds to 
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a regular charge density distribution of flavor branes and we find a family of embeddings of 
D5-branes which reproduces the result of the macroscopic formalism. Finally, in section [5] 
we summarize our results and present our conclusions. 



2 A/" = (2, 2) 2d SQCD from wrapped D3-branes 

The background dual to the two-dimensional gauge theory which we will analyze in this 
section is generated by a stack of D3-branes wrapping a two-cycle of a Calabi-Yau cone 
of complex dimension three, according to the setup: 

CY3 







52 






D3 

























where S"^ represents the directions of the two-cycle C2 that the branes are wrapping and 
A^4 are the directions of the normal bundle to C2. In the above setup a circle represents 
a wrapped direction, whereas the symbols "-" and "." denote unwrapped worldvolume and 
transverse directions respectively. We shall parametrize the cycle C2 by means of two angular 
coordinates [6, (p) with < 6' < vr and < < 27r. The normal bundle A^4 is non- 
compact and will be parametrized by a radial coordinate a, together with three other angular 
coordinates. Moreover, the transverse space will be described by polar coordinates (p. A) 
with < p < 00 and < A < 2tt. The ansatz for the ten-dimensional metric of type IIB 
supergravity corresponding to this setup is: 



H- 



+ 



[w 



+ COS e d(pf ] ) + {dpf + p2 {d\f 



(2.1) 



where dx\ ^ denotes the Minkowski metric in 1 + 1 dimensions, the 
left-invariant one-forms satisfying dw^ 
mass which, for convenience we will take as 



w 



1,2,3) are SU(2) 



i ^^3k /\ rn is a constant with units of 



^ = ^/4:7igs Nc a' 



(2.2) 



with Qs and a' being respectively the string coupling constant and the Regge slope of su- 
perstring theory. The metric ansatz (12. ip contains two functions, z, which controls the size 
of the cycle C2, and the warp factor H. Both of them should be considered as functions 
of the two radial coordinates p and a, i.e. z = z{p,a), H = H{p,a). Moreover, as any 
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other background generated by D3-branes, the ansatz should be endowed with a self-dual 
RR five-form F5. Let us write F5 as: 

F5 = + *J-5 , (2.3) 

where can be represented in terms of a four-form potential C4 as: 

J's = . (2.4) 

The explicit ansatz that we will adopt for C4 is: 

C4 = gxw^ ^'up■ ^ {vf" + cos Bd(\)) AdX + g2^2 /\ {w^ + cos ed(f)) A dX , (2.5) 

with Qi and g2 being two functions depending on both radial coordinates p and a and Q2 
being the volume form of the two-sphere: 

n2 = smOdOAd^ . (2.6) 

We will determine the functions z, H, gi and g2 of our ansatz by imposing the condition 
that our background preserves four supersymmetries. In order to specify this condition it is 
convenient to fix the following vielbein basis of one- forms for the metric (12.11) : 

eO'i = da;°'^ , = ^dO , = ^ sinO dcf) , 

m m 

1 1 



e = — — da , e = — - aw , e = — 5- a w 

zi 2zi 2zi 
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e'^ = ^ a{w^ + cose d(f)) , = H^i dp , = pdX . (2.7) 

2zi 

We shall now impose the vanishing of the supersymmetry variations of the dilatino and 
gravitino of the the type IIB theory. The corresponding Killing spinors e are required to 
satisfy the following set of projections: 

r2356 e = e , r5647 e = -e , roi23(^T-2)e = e , (2.8) 

where e is a doublet of Majorana-Weyl spinors of fixed ten-dimensional chirality and T2 is 
the second Pauli matrix, which acts on the doublet e. In (12. 8p Ta^ a2 ■■■ are antisymmetrized 
products of constant Dirac matrices in the vielbein basis (12.71) . Moreover, the four unknown 
functions z, H, gi and g2 must satisfy the following system of first-order differential equations: 

z , , w?a 



^ = 8"^^ -^{92 - 9i) , z = H , 



pa^ yJZ 

9i=o^H-——H, 91 = -o H + m —{g2 - gi 



z 16 z'^ 



z^ 



g'^ = -]P^H, g, = l^g,-g,)^ (2.9) 

A y/z a 



where the prime (dot) denotes the partial derivative with respect to p (cr). Similarly to other 
backgrounds studied in [23l HEl [31], |l7j, the system (12 .Qp can be reduced to the following 
PDE for the function z(p, cr): 

2pz^{az + z) = (t{pz'^ - 2zz' - 2pzz") . (2.10) 

One can verify that, if z{p, a) is known, the other three functions H, gi and g2 of our ansatz 
can be determined from the BPS system (12.91) . Moreover, if (12. 9p holds, the second-order 
equations of motion of type IIB supergravity are satisfied. Notice also that the Killing 
spinors e satisfy the three independent projections (12.80 . which means that the background 
is |-supersymmetric and four supersymmetries are unbroken. It is easy to check that there 
are two supercharges of each two-dimensional chirality, as it corresponds to the gravity dual 
of an A/" = (2, 2) gauge theory. 

The BPS system (12. 9p can be written in a very compact and suggestive form in terms of 
the so-called (generalized) calibration form /C. For the case at hand /C is a four-form, which 
can be represented in the vielbein basis (12. 7p as: 

/C = j^ Ka,...a,e'''-''' , (2.11) 

where e'^^ ■"'i = e"^ A - ■ ■ Ae"*. The different components of /C are obtained from the fermionic 
bilinears: 

lCa,...a, = H-^ eUT2Ta,...a,e , (2.12) 

where e is a Killing spinor of the background. Taking into account the projections satisfied 
by the spinor and the normalization condition Hi e^e = 1, we get the following expression 
for /C: 

K = A{e''-e''-e^') . (2.13) 

The form K, can be rewritten more compactly if one takes into account that the Kahler form 
of the internal manifold in the frame (12. 7p is just: 

J = e'' - e'^ - e^' . (2.14) 

Then, if Vol(Mini^i) is the volume form of the Minkowski part of the metric (Vol(Mini^i) = 
one can write the calibration form /C as: 

JC = Vol(Mini,i) A J . (2.15) 

One can verify that the supersymmetry preserving conditions (12. 9 p can be rewritten as: 

dlC = *dC4, d(^H-^ * IC^ = , (2.16) 

where the star * denotes the Hodge dual in the ten-dimensional geometry (12. ip . Notice that 
the first of the conditions (12.160 implies the following relation between the RR five-form field 
strength and the calibrating form: 

= die + * die. (2.17) 
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2.1 Integration of the BPS system 



The brane setup studied in this section can be reahzed in the context of five-dimensional 
gauged supergravity [5]. In this approach one formulates an ansatz for the 5d metric, scalar 
and gauge fields in which the different functions depend on one radial variable and, therefore, 
the corresponding BPS equations are ordinary differential equations. This BPS system was 
obtained in section 3.2 of [5] and will not be repeated here. After uplifting the ansatz of 
[5] to ten dimensions, and by performing a suitable change of variables, one can define two 
radial coordinates p and a in such a way that a metric of the type (12. ip and a five-form as 
in (I2.3l) -( l231) are obtained (see refs. [Ml Il7] for similar analysis in other brane setups). 

Although we have not been able to analytically integrate the BPS system of [5], it turns 
out that the different first-order equations can be combined to produce a single second-order 
differential equation, whose solutions allow one to solve our BPS system (12. 9p in implicit 
form. In order to present this solution, let us introduce an auxiliary function r = t{z), 
defined as the solution of the following ordinary second-order differential equation: 



2e'^ - — 



dz 



(2.18) 



Then, one can show that the PDE (I2.10p can be solved in implicit form by taking z{p, a) as 
the solution of the equation: 



o2t{z) 



+ 



a 



dr 

dz 



(2.19) 



The other functions of our ansatz can also be obtained. Indeed, the warp factor H{p, a) 
corresponding to the function z{p, cr) of (I2.19P is given by: 



H 



ze' 



-2t 



(2.20) 



while the functions gi{p,a) and g2{p,<j) are: 



a 



92 



8 ^/z \dz/ 



(2.21) 



One can verify that z, H, gi and g2, as given in eqs. (I2.18p -( 12^2T|) . satisfy the system of 
PDEs written in (12. 9p . It is also easy to check that, when m is given by (12. 2p . the functions 
(I2.2ip give rise to the correct quantization of the flux of for a stack of Nc D3-branes. The 
auxiliary equation (I2.18P can be easily integrated numerically. This numerical solution can 
be used in (I2.19p -f l2.2ip to get the different functions of our ansatz in the (p, cr) plane. Before 
performing this numerical analysis, let us present in the next subsection some approximate 
results for the region of the background in which the value of the function r is large which, 
as we shall argue, will turn out to correspond to the UV region of the corresponding dual 
field theory. 
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2.1.1 Approximate UV solution 

For large r the master equation fl2.18p can be written as: 

d'^T o-r- dr d 



which can be immediately integrated, namely: 

^ = e'^ + c, (2.23) 

dz 

with c being a constant of integration. An additional integration gives r as a function of z: 

= e^.(.,-.) _ 1 ' (2-24) 

where z^, is a new constant of integration which represents the value of z for which r — > cxd. 
By making use of fl2.24p in f l2.19p . we arrive at the following implicit relation for z{p,a) in 
the large r region: 



P 



2 



2 



(2.25) 



Notice that for large r one has z ~ z=k and, according to (12.251) . we are far from the origin in 
the (p, a) plane. Thus, large r corresponds to the UV region of the background, as stated 
above. Furthermore, in this deep UV region we can take e^"^'*-'^ - 1 ^ 2c(z* — z) on the 
right hand side of (12.251) and, therefore z{p, a) can be approximately written as: 



z^z.~fi . (2.26) 

Moreover, it follows from (I2.19P that in the deep UV where, at leading order, z ^ z^, and 
^ f» e^"^, one has: 



a' + y^p^ ^ ^e'^ . (2.27) 



m 



On the other hand, using these UV estimates in (I2.20p and (I2.2ip . one gets: 
^ ^ 1 _ _l 1 



m4 ( a2 + p2 )2 ' 8^4(^2 + ^^2)2' o'^ \ p^ ' 

(2.28) 

The above results suggest that, in the UV, the different quantities depend on o and p through 
the combination + ^Jzl (? . Accordingly, let us introduce the following new variables u 
and a: 

z^)^ p = u cos a , (2.29) 
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with < M < oo and < « < 7r/2. Since the combination cr^ + ^Jzl is just m^, we get 
that: 

^ (2.30) 



Using these resuhs for H and z, one immediately verifies that the UV metric becomes: 



ds 



m 



uv 



u 



m 



+ 



1 (du) 



+ 



+ - 



{daf + 



sin^ a 



w'Y + {wy + («;•' + COS ^ci0)' + cos'a(ciA)' . (2.31) 



Notice that the first fine corresponds to the metric of an AdS^ space, with two of its directions 
compactified on an S"^. The second fine is the metric of a five-sphere fibered over the S"^ . 



2.2 Probe analysis 

In order to make contact with the gauge theory dual to the background just studied, let 
us consider a probe color D3-brane. This probe brane is extended along (x°,x^,^,0) at 
fixed values of the other coordinates. We will first assume that the worldvolume gauge fields 
of the D3-brane are not excited. By explicitly computing the DBI and WZ terms for this 
configuration by using the metric (12. ip and the RR four-form potential corresponding to the 
F5 of (I2.3p -( l231) (F5 = dC^), one realizes that the D3-brane is only at equilibrium if the 
coordinate a vanishes. Moreover, one can check that a = is the SUSY locus of the color 
D3-branes. The calculations leading to these conclusions are the same as those performed 
in section 4.1 of [16] for the case with eight supersymmetries and, therefore, they will not be 
repeated here. 

Let us next consider a color brane located at the no-force point a = and let us assume 
that we switch on the worldvolume gauge field F in such a way that its only non-vanishing 
components are those along the unwrapped directions x^. By expanding the DBI action in 
powers of F and looking at the coefficient of the F"^ term, we get the holographic expression 
of the Yang- Mills coupling Qym- Again, the calculation is the same as in [16] and gives the 
result: 

1 z{p,a = Q) 



9ym gsm"^ 



(2.32) 



Notice that the right-hand side of (I2.32p is proportional to the size of the cycle, which is 
controlled by the function z. Moreover, since z{p,a) in fl2.32p is evaluated at a = {i.e. at 
the SUSY locus), the remaining radial coordinate p plays the role of a holographic coordinate. 
In the next subsection we will present numerical results for z{p,a). It is however possible 
to extract relevant physical information by using the analytical results of subsection 12 . 1 . 1 1 in 
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Figure 1: On the left we plot z, obtained from fl2.18p and f l2.19p . as a function of the 
dimensionless variables mp and ma. We have fixed the UV constants and c of fl2.25p to 
the values 2* = 2 and c = — 1/2. On the right we represent the warp factor H for these same 
values. We only show the region where H is monotonic. 



the UV region of the geometry. Indeed, by using fl2.26p for a = and large p we get: 



9\m 9i '2m^9sp'^ ' 

with gl = rn^gs/z^.. It follows from fl2.33p that gyM approaches the fixed-point value gyu = 
g^ in the UV and deviates from it by means of a power law in the holographic coordinate p. 
Notice that the running with a power of the right-hand side of fl2.33p was to be expected 
from dimensional arguments. It is just the behavior obtained in one-loop perturbation theory 
if p is taken to be proportional to the renormalization energy scale p. The minus sign of the 
running term in f l2.33p is also expected: it just means that gyM decreases when we move 
towards the UV. However, the coefficient of this term does not coincide with the one found in 
perturbation theory if the naive identification p = 2Txa' p is used. This fact is not surprising, 
given the amount of SUSY preserved by our background. 



2.3 Numerical results 

The gauged supergravity solution fl2.18p - fl2.2ip of the BPS system (12. 9p can be easily evalu- 
ated numerically. One must first solve the master equation fl2.18p and use the function t{z) 
obtained in this way to find z(p, cr) by solving the implicit equation f l2.19p . The warp factor 
H and the functions gi and (72 are then straightforwardly obtained from fl2.20p and fl2.2ip . 

The result of this numerical analysis for ^(p, cr) and (p, a) has been plotted in figure [1] 
We notice from this figure that z(p, a) grows monotonically when we move away from the 
origin of the (p, a) plane and reaches a constant asymptotic value when p^ + a"^ — )■ 00. On the 
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contrary H{p^ a) increases as + o"^ decreases until it reaches a maximum and then starts 
to decrease. This non-monotonic behavior of the warp factor is opposite to the one expected 
for a good holographic dual of a gauge theory and it indicates to us that the supergravity 
solution is not trustable in this IR region. Notice that this problem in the IR is shared by 
other backgrounds constructed from branes wrapping cycles. 

2.4 Addition of flavor 

In this section we will modify the previous setup to include the effect of degrees of freedom 
corresponding to fields in the fundamental representation of the gauge group (flavors). This 
can be achieved by including in the setup a new set of D3-branes which introduce a new 
open string sector. These flavor D3-branes should be extended along the two Minkowski 
directions, as well as along a two-dimensional non-compact submanifold of the normal bundle 
A^4. Moreover, the flavor branes should be located at a fixed point of the two-cycle C2 and 
of the transverse R^. The corresponding brane array is: 

CY3 
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In order to determine the particular embedding of the flavor D3-branes in the normal bundle 
A^4 we will require that the configuration described by the array written above preserves the 
same four supersymmetries as the unflavored background. In section [2.4. II we will introduce 
a general family of such embeddings, which can be briefly described as two-planes embedded 
in A^4 (which has the structure of R^ fibered over the two-cycle C2). Notice also that the fixed 
value pQ of the p coordinate of the flavor branes represents the distance between the two 
sets of branes in the transverse R^, which should be related to the mass mq of the matter 
fields as mq ~ Pq/a' ■ 

When Nf « N^ the flavor branes can be treated as probes which do not alter the 
unflavored supergravity background. This is the so-called quenched approximation. Here we 
will concentrate in the limit in which Nf is large and of the same order as N^ and, therefore, 
the effects of the backreaction must be taken into account. Indeed, when Nf ~ A^^^ one must 
deal with a coupled gravity plus branes system in which the branes are dynamical objects 
that act as sources for the different supergravity fields. In particular, the WZ term of the 
D3-brane worldvolume action couples to the RR four-form potential C4 in the form: 

Swz = T3J2 [ ^4 , (2.34) 

where the hat over C4 denotes its pullback to the D3-brane worldvolume Ai^. In general 
a term like fl2.34p is a source for the RR five-form F5 which gives rise to a violation of 
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its Bianchi identity on the flavor brane worldvolume. Actually, solving the corresponding 
equations of motion of the coupled system when the flavor branes are embedded along a 
given fixed submanifold Ai^ is a formidable task which, in practice, is not possible to tackle. 
For this reason we will follow the approach of ref. |22] (see [31] for a review) and we shall 
distribute the A^/ — )• oo flavor branes along a continuous set of mutually supersymmetric 
embeddings. In this approach one performs the following substitution in the WZ term: 

J2 [ ^ f ^] A C4 , (2.35) 

Jm4 J Mio 

where is a six-form (the smearing form) which encodes the RR charge density associated 
to the continuous distribution of flavor branes. Actually, is just the source in the modified 
Bianchi identity of F5, namely: 

rfFs = 2k?oT3(] , (2.36) 

where 2k^q = (27r)^ (7^ {(^'Y- It is clear from (12.361) that, in order to include the backreaction 
effects, we must modify our ansatz fl2.3l) - fl2.5p for F5. Actually, instead of (12. 4p we shall take 
J-5 as given by: 

^5 = rfC4 + /5 , (2.37) 

where C4 is still given by the ansatz (12.51) and /s is a five-form which incorporates the 
modification of the Bianchi identity (I2.36p . namely: 

dh = 2t^l^T^n. (2.38) 

As our flavor branes will be located at a fixed value of p = pg, it follows that the smearing 
form VL should contain a S{p — pq) in its expression. Moreover, it is clear from (12.381) that fl 
is a closed form, i.e. dQ = 0. These two conditions are satisfied if Q is of the form: 

2kIq T3 = 6{p - pq) dp AdA , (2.39) 

where A is a four-form depending on a and on the angular coordinates. Eq. (I2.39P implies 
that /s can be taken as: 

/s = Qip-pQ)dA, (2.40) 
and, therefore, the total J-5 is given by: 

= dC^ + Q{p-pQ)dA . (2.41) 

The precise form of A (and therefore of Q) is obtained by computing the charge density that 
results after averaging over a particular family of equivalent embeddings that are mutually 
supersymmetric. Instead of following this "microscopic" procedure to determine A, we will 
follow the "macroscopic" approach of [SOj [17] and we shall try to get the form of A by 
requiring the preservation of supersymmetry and the compatibility of the charge density 
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with the metric ansatz (12. ip . Later on, in subsection I2.4.H we will precisely characterize 
the particular set of embeddings that give rise to the charge density determined in the 
macroscopic formalism. 

In the flavored case we shall adopt the same ansatz fl2.5p for the RR potential C4 as in the 
unflavored background. Moreover, from the expression fl2.4ip of J-5 the similarity between 
C4 and A is quite obvious. For this reason it is quite natural to adopt an ansatz for A with 
the same structure as in (12.51) . namely: 

A = Li(o-) Aw"^ A {w'^ + cos 6*6/0) AdX + L2{cr) Q2 A {w^ + cos Odcj)) A d\ , (2.42) 

where Li and L2 are functions of the coordinate a to be determined. Since dA is given by: 

dA = {L2-L1) n2Aw^ Aw^ AdX + (^Liw^ Aw^ + L2 ^^2 ) A c/d A (w^ + cos 6*^0 ) A rfA , 

(2.43) 

it follows that the flavored BPS ansatz for J-'s can be obtained from the unflavored one after 
performing the following substitutions: 

gi-^ gi + Li e(p - pq) , (2 = 1,2), 



g2- gi-> g2- gi + {L2 - L^) e(p - pq) . (2.44) 

One can now repeat the supersymmetry analysis for the new ansatz. Obviously the flavored 
BPS system is obtained by performing the substitutions (I2.44p in (12. 9p . namely: 



z' = Sm^ 



H 



[g2 - gi + {L2 - Li) e(p - Pq) ] 



z 

8 16 ^2 

gi = -I— H' + m^^{g2 - gi + {L2 - U) Q{p - pq)) - L^Q{p-pQ) , 

Z Z2 

1 pa^ 

g-i = — , 

g2 = -{g2 - gi) - [L2 - -(L2-Li)]0(p-pq) . (2.45) 
cr o 

By analyzing the compatibility of the different equations in the system (I2.45P one discovers 
that, in general, the (p, ex) crossed derivatives of the functions g\ and g2 are not equal. Indeed, 
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one can prove that: 



dp9i - d^g'i = - dpg2 = --[2(^2 - Li) - aL2]6{p - pq) . 



(2.46) 



Therefore, in order to get rid of this unwanted singularity we must require the following 
condition on Li and L2: 

-{L2-L1). (2.47) 



a 



Notice that, after imposing fl2.47p . the last equation in fl2.45p is the same as the one in the 
unflavored system (12. 9p . Furthermore, one can verify that any solution of the BPS system 
(I2.45P and the compatibility condition (I2.47P also solves the equations of motion of gravity 
with source terms coming from the smeared branes. Moreover, in this case one can also write 
a single PDE for z{p, a), which now has a source term parametrized by L2 — Li, namely: 



2pzy/z{az + z) = a {pz'"^ - 2zz' - 2pz z") + 



16m^(L2 — Li) Z2 



a 



6{p-pq). (2.48) 



The compatibility condition (I2.47P allows one to determine ^2(0') in terms of Li{a). In 
principle, there are and infinite number of solutions to (12.470 . However, it turns out that 
there is a particularly simple solution which satisfies certain regularity conditions which 
should be satisfied on general physical grounds. Since the complete analysis is very similar 
to the one performed in section 3.3 of [ITj, it will not be repeated here. Instead, we will write 
down the solution and we will verify that it satisfies some of these physical requirements. 
The solution of (I2.47P we will be interested in from now on is the one obtained by taking 
Li((t) = 0, namely: 

Li = , L2 = Ca^ , (2.49) 

where C is a constant. By using these values of Li and L2 in eqs. (12.431) and (I2.39P one can 
readily get the corresponding expression of the smearing form Q, i.e.: 



n 



C 



2kI,T, 



6{p - Pq) dpAdXAQ2 A 



cr^ Aw'^ + 2a da A 



(2.50) 



which, in terms of the vielbein basis (12.71) . takes the following form: 

Am'C 6{p - pq) 



n 



2«?o^3 pVni 



(2.51) 



An important test that the smearing form Q must pass is the regularity and positivity of the 
mass density distribution of the flavor brane sources. This mass distribution density can be 
obtained by looking at the smeared version of the DBI action of the flavor branes. Indeed, 
taking into account that, for a SUSY conflguration, the induced volume form is just the 
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pullback of the calibration form /C and that the smearing is performed by taking the wedge 
product with Q, this DBI action takes the form: 

g flavor = / A /C . (2.52) 

JMio 

Thus, QAK. can be interpreted as the mass distribution of the system of flavor branes which, 
being a ten-form in a ten-dimensional space, is proportional to the volume form Vol (A^io) 
of the ten-dimensional manifold. Actually, from f l2.5ip and the expression of K, in (12.131) one 
straightforwardly gets: 

The function multiplying Vol (A^io) in (I2.53P represents the mass density of the p = pq slice 
of the space in which we are smearing the flavor branes. The positivity of this mass density 
implies that the constant C must be negative. In the next subsection we will flnd a family 
of supersymmetric embeddings which gives rise to the charge distribution (I2.50p and we will 
identify the constant C with the parameters of this set of smeared branes. 

2.4.1 A microscopic interpretation 

In order to characterize a family of embeddings of the flavor D3-branes that preserve the 
supersymmetries of the background, let us parametrize the one-forms in terms of three 
angles 9, and if) as follows: 

= cos 4' d9 + sin sin 9 d(l> , 
w'^ = sin ip d9 — cos ip sin 9 dip , 

w'^ = dtp + cos 9 d(j) , (2.54) 

with 0<^<7r, O<0< 2n, < < An. Let us next introduce the coordinates 
(z = 1, ■ ■ ■ ,4) as: 

+ ty^ = a COS (^) ('^+^) , y'^ + ty^ = a sin (^) (^-^) . (2.55) 
One can verify by direct calculation that: 

J2{dy^Y = {dar + ^J2{w^)\ (2.56) 

i i 

which shows that the normal bundle A4 has the structure of flbered over the {9, (p) two- 
sphere and that the ?/*'s are just the cartesian coordinates of this four-plane. Moreover, one 
can also check that: 

dy^ A dy^ + dy^ A dy'^ = Aw"^ + 2a da A . (2.57) 
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Using fl2.57p in f l2.50p one immediately proves that the smearing form Q can be written in 
terms of the coordinates as: 

2C 

n = 2— 6{p - pq) rfp a dA a ^2 a [ rfy^ A dy"^ + dy^ ^dy'^] . (2.58) 

The standard way to determine if a given D-brane embedding is supersymmetric is by verify- 
ing the fulfillment of the kappa symmetry condition on the Killing spinors of the background. 
Equivalently, one can demonstrate the supersymmetric nature of the embedding by checking 
that the puUback of the generalized calibration form on the D-brane worldvolume A^4 is 
equal to the induced volume element. In our case this calibration condition reads: 

t = Vol(A^4) , (2.59) 

where /C is the four-form written in f l2.13p . According to the brane array in the beginning 
of subsection 12.41 let us consider a D3-brane which sits at a particular point of the {6, 0) 
two-sphere and that is localized in the plane parametrized by (p, A). In addition, the D3- 
brane is extended along a codimension two surface of the four-dimensional plane spanned by 
the y^ coordinates. The puUback of /C for these embeddings takes the form: 

/C = -^dx° A dx^ A [dy^ A dy^ + dy^ A dy^] , (2.60) 



where the pullback of the coordinates is understood. Moreover, let us assume that the 
embedding in the normal bundle A^4 is described by linear relations of the type: 

y^ = a,y^ + h, y^ = a2y^ + b2, (2.61) 

which represent a two-plane in M^. For these linear embeddings one has: 

^ = -/^ ( 1 + ai 02 ) c/x° A c/x^ A dy^ A dy'^ . (2.62) 



On the other hand, the induced metric on the D3-brane worldvolume can be written as: 

ds^ = H-^'^dxl^ + [(1 + a?) {dy^f + (l + a^) [dy^f] , (2.63) 

^/ Z 



and the corresponding volume form is 



Vol(7W4) = ^\ {l + al){l + al)dx^ ^dx^ ^dy^ ^ dy^ . (2.64) 



It is now straightforward to prove that the calibration condition holds if Oi = 02. Let us 
parametrize: 

0-1 = 02 = —cot 7 , (2.65) 
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where 7 is a constant. Then, the embedding in the ?/* hyperplane is characterized by the 
equations: 



fi = cos'jy^ + siwjy^ - ci = , 

/2 = cos 7?/^ + sin-f - C2 = , (2.66) 

where ci and C2 are new constants. Eq. (12.661) defines a family of embeddings parametrized 
by three parameters (7, ci and C2)lll Notice that changing 7 by 7 + vr is equivalent to taking 
Ci — —Ci in (12.661) . Thus, we will take 7 in the interval < 7 < vr. 

Let us now compute the charge density six-form Q that results after averaging over this 
family of embeddings. We will homogeneously distribute the D3-brane embeddings over the 
{6, 0) two-sphere and we will consider branes with a given value pq of the radial coordinate 
of the transverse and distributed over the angular direction parametrized by A. Thus, Q 
can be written as: 

^ r / X , (iA sin 9d9 A dd) ^ . „x 

n = 5{p-pQ)dp A — A ^ AT , (2.67) 

where we have normalized the distributions in the coordinates A and (6', 0) by dividing by 
the volume of 5*^ and 5*^ respectively. In (12.67^ F is the two-form density obtained after 
averaging over the family of planes (I2.66p . It can be written as the integral of the volume 
element of the complement of the two-planes (I2.66P in iV4, namely [28]: 



j dp6{h)6{f2) [dhAdh] , (2.68) 



where dp is an integration measure depending on the parameters 7, ci and C2 and the exterior 
derivative only acts on the y^ variables of the fj. The natural integration measure dp for the 
family of planes (I2.66P is: 

d^ 

dp = Uf — dcidc2 , (2.69) 

TT 

with rif being the density of fiavor branes. The integrals over ci and C2 in fl2.68p can be 
immediately performed by making use of the two delta functions . Moreover, since: 

dfi A df2 = cos^ 7 dy^ A dy^ + sin^ 7 dy^ A dy^ + sin 7 cos 7 [ dy^ A dy'^ - dy'^ Ady^] , (2.70) 

we get: 

r rf7 [dfi A ci/2] = ^ [dy' A dy^ + dy' Ady'] . (2.71) 
Jo ^ 

^Actually, there is a much larger family of calibrated embeddings for this background. If we complexify 
the coordinates as 2;^'^ = y^'^ + iy^'^, then any submanifold in N4 defined by a holomorphic relation of 
the type — f{z^) satisfies (|2.59p . In particular, a linear relation such as az^ + f3z'^ = constant defines 
a complex line in which generalizes (|2.66|) . However, this more general family of planes gives rise to the 
same fl as the one obtained from p.66|) and, thus, we will not consider any of these more general embeddings. 
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Figure 2: Plots of z{p, a) for the flavored background of D3-branes. The flavor branes are 
located at m pq = 1.8 for two different values of x = . The plot on the left (right) 
corresponds to x = 2 (x = 1/2). The solution for p < pq is the unflavored one with UV 
parameters 2;* = 3 and c = —1.15. 



Thus, the two-form T can be written as: 

T = ^ [dy'Ady' + dy' A dy'] , (2.72) 

and the resulting Q is, indeed, of the form f l2.58p with the constant C related to the density 
Hf as: 



2.4.2 Numerical results 



Let us now present the results obtained by the numerical integration of the BPS system 
([235]). As in [231 ini EH IT], our strategy will be to integrate the PDE equation fICTD 
for z{p,a). Actually, after using the values of Li and L2 written in f l2.49p and the relation 
between C and n/ displayed in eq. ( ]2.73p . the PDE ( 12.48^ takes the form: 

2pzy/z{az + z) = a {pz'^ - 2z z' - 2pzz") - — ^ o z^^ 5{p - pq) . (2.74) 

Since the system f l2.45p reduces to the unflavored one for p < pq, we can assume that 
the solution of (I2.74p is given by f l2.18p -( l2.19p in this region. At p = pq the function z is 
continuous while, according to the first equation in f l2.45p . z' has a discontinuity given by: 



z{pq + e,(T) - z{pq-e,a) = - \ ^ . 2.75 

2m^ Nc pq 

We have integrated f l2.74p by using the continuity at p = pq and the change in derivative 
f l2.75p as boundary conditions. The results are shown in figure |2]for two values of Uf. In our 
numerical analysis we have been able to integrate 2; up to a finite value of p, after which the 
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solution becomes highly oscillatory and unstable. As is evident from figure [21 the function 
z has a wedge shape at the position p = pg of the fiavor branes. Thus, there will be a 
curvature singularity at this point, which is actually required by the Einstein equation in 
order to match a similar term coming from the energy-momentum tensor of the fiavor brane 
source. In general, the addition of fiavor makes z(p, a) grow slower or to decrease as we move 
towards the UV for p > pg. In view of the relation fl2.32p between z and the YM coupling, 
this behavior was to be expected since, contrary to what happens with the non-abelian color 
fields, the fiavor matter fields always make the coupling grow in the UV. Thus, our fiavored 
geometry correctly encodes this non-trivial property of the renormalization of gauge theories. 



3 A/" = 2 3d SQCD from wrapped D4-branes 

In this section we will analyze some gravity duals of three-dimensional TV = 2 supersymmetric 
gauge theories which are generated by D4-branes wrapped along two-cycles of a Calabi-Yau 
threefold. Our brane setup, including fiavor, is very similar to the one studied in section |2l 
namely: 

CY3 

. ^ 
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We will be working in the context of type IIA supergravity and we will adopt the following 
ansatz for the metric and dilaton in the string frame: 

z 

IV? 

^2 



dxlo + A,((de? + sm^e(d(p?) 



+ 



\ [ {daf + y [ (^1)2 + + + cos ^ ] j + (dp)' 



Z2 



where —00 < p < +00 and the parameter m is given by: 

^ = {87ig,N,ya' . 
m 



(3.1] 



(3.2) 



The similarity of the metric in (13. ip and the one corresponding to the D3-brane case in (12. ip 
is manifest. However, notice that, in this case, the background should be endowed with an 
RR four- form field strength F4. In analogy with the D3-brane case, we will take F4 to be of 
the form: 

Fi = dCs + Q{p-pQ)dA, (3.3) 
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with C3 and A being given by: 



C3 = giw-^ Aw^ A {w^ + cos 9d(f)) + g2^2A [w^ + cos 

A = Li((t) Aw"^ a {w^ + cos Odcj)) + ^2(0-) n2 A {w'^ + cos 6*6/0) , (3.4) 

where fl2 is the volume form of the 5*^ written in fl2.6p and gi{p, a) and Li{a) are functions 
to be determined. Notice that we have aheady included in (13. 3 p the term that modifies the 
Bianchi identity for F4 due to the presence of the flavor D4-branes at p = pg. Indeed, by a 
simple calculation from (13. 3 p we get: 

dF4 = 5{p - pq) dp AdA . (3.5) 

The source term on the right-hand side of (13. 5p comes from the coupling of the D4-brane to 
the RR fields through the WZ term of the worldvolume action, Swz = T4 J Q A C5, where 
fl is the RR five-form charge density of the flavor brane distribution and C5 is the five-form 
potential of Fg ( F4 = * Fq). The corresponding Bianchi identity with source is: 

dFi = 2KloTin , (3.6) 

which allows one to identify Q for our ansatz with: 

2kIq T4 n = 5{p - Pq) dp AdA . (3.7) 

In order to find the system of BPS equations for a type IIA background with metric, dilaton 
and RR four-form as in our ansatz, let us consider the following vielbein basis for the metric 



m ' m 



c iJ4 Hi 1 7 H* 2 

e = — r dcr , e = — -aw , e = — -aw , (3.8) 

zi 2zi 2zi 

}jl 1 
= — J a [w^ + cos 9 dcf) , = dp . 

2zi 

We will now impose the following projections on the Killing spinors of type IIA supergravity: 

r3467 e = r3458 e = e , T 56789 e = (3.9) 

where the indices of the constant Dirac matrices refer to the basis (13. 8p . Notice that the 
projections (13. 9p leave four supercharges unbroken, corresponding to a three-dimensional 
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theory with Af = 2 supersymmetry. The corresponding system of BPS equations is: 



) m 



m? a 



92-91 + {L2 - Li) e(p - 



H. 



91 



9'2 



92 



'z 
1 



H 



m a 



1 ^3 H 
— - — H + m 



» z 
4 



3 

Z2 



92 - 91 + {L2 - Li) Q{p - pq) -LiQ{p-pq) 



H. 



2 . 

- [92 

a 



9i) 



L2 — {L2 - Li) e(p - pq) . 



a 



(3.10) 



Combining the equations of fl3.10p in different ways we can check that the cross derivatives 
of both 9i and 92 depend on the ordering of the differentiation. In order to avoid this 
unwanted singularity we need to impose again fl2.47p . which imphes that the Lj's drop from 
the last equation in f lS.lOp which, therefore, is the same in the flavored and unflavored cases. 
Moreover, any solution of f l3.10p with Li and L2 satisfying f l2.47p also solves the equations 
of motion of the gravity plus (smeared) branes system while, as in the D3-brane case, we 
can write a single PDE for the function z{p,a), which has a source term parametrized by 
L2 — Li, namely: 

16m^2;^/^ 



2 z \fz [az + i) = o" ^z'^ — 2z z" 



a 



{L2-L1) 6{p- Pq). 



(3.11) 



As in the D3-brane case, the BPS system fl3.10p can be recast in terms of the appropriate 
generalized calibration form /C. In this case /C is a five-form, defined in terms of bilinears of 
the Killing spinors as follows: 



(3.12) 



where e is normalized as if 4 e = 1 . After using (13. 9p we get the following expression for 
/C: 



A{e 



,58 



„67 



(3.13) 



Using the expression f l3.13p for /C, as well as the ansatz fl3.3p -( l3^ for F4, we can easily verify 
that the following two conditions: 



d 



IC 







d mic 



(3.14) 
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reduce to the BPS system f l3.10p . In (13.141) the Hodge dual is taken with respect to the string 
frame metric in (13.11) . 



3.1 The unflavored solution 



In the unflavored case, the setup we are studying here can be reduced to the one analyzed 
in [5], where a configuration of M5-branes wrapping a two-cycle and preserving four super- 
charges was considered. By compactifying one of the worldvolume directions of the M5-brane 
and dimensionally reducing along it, the setup of [5] reduces to the one studied here. In [5] 
seven-dimensional gauged supergravity was employed to find a system of first-order BPS 
equations for the functions of the ansatz. These functions depend on a single radial variable 
and, therefore, the BPS system of [5j is a system of ordinary differential equations. As in 
section [2], one can perform a change of variables which converts the ansatz of [5] (uplifted to 
ten-dimensions) into our two-variable ansatz. Moreover, this relation can be used to find a 
solution of the BPS system (I3.10p with Li = L2 = from a solution of the BPS system in 
[5]. Indeed, let t{z) be the function that solves the following differential equation: 



— 2^— V 

dz"^ \dz J 



2 

dz 



4e2- - — 



1 

Xz 



(3.15) 



It turns out that we can write a solution of the unfiavored BPS equations in terms of this 
function t{z). First of all, one can verify that the PDE (13. lip with Li = L2 = is solved 
by a z(p, a) that is given implicitly by means of the following equation: 



e^^ ^ \ ~ Am? ' 

Moreover, the warp factor H is given in terms of z(p, a) and t{z) by: 



(3.16) 



H 



4a2e2- + p2yi(g)%-2r 



(3.17) 



while the functions gi and g2 are: 



9i 



pH /drV 



QAm? \dz J 



d7 



dz 



92 



pe 



-2t 



(3.18) 



One can easily show that the RR fiux for the solution (13.180 is correctly quantized when the 
constant m is given by (13. 2p . 

It is now straightforward to solve numerically first (I3.15P for t{z) and use this result in 
the implicit equation (I3.16p to get z{p, a). The results are plotted in figure [3] for z(p, a) and 
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Figure 3: Plots of z and H for the unflavored background of D4-branes. We are taking the 
UV parameters 2;^, = 2 and c = 1. On the right plot we only show the region in which the 
warp factor H is monotonic. 



H{p,a). They are very similar to the ones found in section [2] for the wrapped D3-branes 
and all the comments made in subsection 12.31 also apply to the present case. In particular, 
the function H is not monotonic in the IR and one should only trust the region in which H 
decreases as we move towards the UV. 



3.1.1 The approximate unflavored solution in the UV 

When T — 7- 00, z approaches a constant value and, therefore, the term l/z on the right-hand 
side of fl3.15p becomes negligible compared to e^"^. Thus, approximately, the master equation 
fl3.15p becomes: 

(Pt f dT\^ r. dr , , , , 

— + 2 - ^Se^--, (r^oo). (3.19) 

az'^ \dz J dz 

This equation can be integrated once to give: 

^:^2(e- + ce-), (3.20) 

while an additional integration gives r as a function of z, namely: 

= y^cot [4v^(z, - z)] , (3.21) 

with c and z^, being constants. Notice also that, for large r, we can take at leading order 
^ ~ 2e^^ and the implicit equation fl3.16p becomes: 

2 _ 1 1 
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Eqs. (13.211) and (13.221) show that for large r one approximately has z ^ z^, and we are far 
from the origin in the (p, a) plane. Moreover, by performing an expansion around z^, in 
(I3.22P and by keeping only the first non-trivial term we have: 



z{p,a) 



1/4 



(3.23) 



Notice also that, in the far UV where 2; ^ , one has 



4^2 



(3.24) 



Using these approximate UV expressions in (13.171) we arrive at the following UV estimate of 
the warp factor H: 

' 1 



E 



2m 



3 ' 



(3.25) 



whereas from (I3.18P we get that g\ and are approximately given by: 



9\ 



PZ: 



?>2m? 



p^ z'J'^ + cr^ 



^2 J 



92 



pZ: 



1/4 



32m z^^ + a^ 



(3.26) 

The above analysis suggests that in the UV region the combination p^z'J'^ + 0"^, appearing in 
the denominator of both z and if, plays a significant role. Having this in mind let us define 
a new set of variables u and a as in (I2.29P but now with < a < vr. In general, z and H 
depend on both coordinates, u and a, but in the UV region the a dependence disappears. 
Actually, their expressions when m — )■ 00 are: 



z ^ z^ 



^1/2 



2^2 \^ 



2'^ 



1/4- 



3/2 



mu 



(3.27) 



Using these values in the metric ansatz (13. ip we end up with the following expression: 



2^2 ( 



mu 



3/2 



+ 



\ 1/4 
1 



2^/2 zV^ ^^^^ 



da^ + - sin^ a 
4 



+ 



2V2zV 



(mu 



,3/2 



+ 



COS 



,(3.28) 



where dfl"^ = + sin^ 9 dcj? is the line element of the {6, 0) two-sphere. In order to interpret 
the meaning of the results just found, let us recall that, given a background of type IIA 
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theory such as the one we are considering, one can generate a solution of eleven-dimensional 
supergravity by uplifting the metric in terms of the standard formula: 



dsf, = e-^'^dsio + e^'^idx'Y , (3.29) 

where is the eleventh M-theory coordinate. We shall apply fl3.29p to the ten-dimensional 
UV metric and dilaton written in eqs. 03.281) and (13.271) . After changing the radial variable 
M by a new coordinate y, defined as: 

n 2m 



y 



(3.30) 



the resulting eleven-dimensional UV metric becomes: 



dsli 



dx\ 3 + — ^ dVt^ 



+ 



4m^ 



do? H — sin^ a 
4 



y 



+ 



[w 



1\2 



+ {w'^f + {w^ + COS 



(3.31) 



where dxl^ = dx\2 + {dx^Y- From fl3.3ip we conclude that the uplifted metric is of the 
form AdS-j x S*^, with the AdSj having two of its directions compactified in a two-sphere and 
with the S'^ being fibered over this S"^ . Notice that the radius of the AdSj {S^ ) factor in 
the asymptotic metric (I3.3ip is 1/m (l/(2m)). These results are, of course, consistent with 
the origin of the solution of [5J, as coming from M5-branes wrapping a two-cycle. 



3.2 Probe analysis 

In complete analogy with the analysis performed in section 12.21 we can test our background 
by using a color D4-brane extended along the Minkowski directions and those of the two- 
cycle. Again, the equilibrium position of such a brane occurs at the no-force SUSY locus 
cr = 0. The coupling of the Yang-Mills theory living on the color branes can be obtained 
by introducing a worldvolume gauge field along the unwrapped directions and expanding at 
second order in the gauge field strength. One gets: 

' - "^^'^ = °^ . (3.32) 



From our numerical results of figure [3] we notice that z{p, a) grows monotonically as p 
increases and reaches a constant asymptotic value when p ^ oo. Using this result we 
conclude that qym decreases as we go to the UV. We can have an analytic estimate of this 
behavior by using the asjTiiptotic results of subsection 13.1.11 Indeed, plugging (I3.23P in 
(I3.32P we arrive at: 

1111 , , 

(3.33) 
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where is the asymptotic value of the YM couphng (g^ = lixg^rr?^^' j z^. Again, the 
behavior displayed by our holographic model coincides qualitatively with the one found in 
the perturbation theory of the gauge theory dual. Indeed, if one relates the coordinate p 
with the energy scale of the field theory by means of the naive relation p = 2Tia' p, one gets 
that the right-hand side of (13.331) runs as p~^, which is the same running obtained with the 
one-loop beta function (although the numerical coefficients are different). 



3.3 Flavored solutions 



Let us now analyze the system (I3.10p in the flavored case. According to our macroscopic 
analysis, the source flve-form Vt is parametrized by two functions Li and L2, which must 
satisfy the consistency equation (I2.47p . However, as was the case for the D3-brane system 
of section [2], there is a simple solution to (I2.47p . namely (I2.49p . which leads to physically 
sensible results. Therefore, from now on we will assume that Li and L2 are given by the 
expressions displayed in ( ]2.49p . The corresponding smearing form f2 is: 



C 



5{p - pq) dpAn2 A 



a^w^ Aw^ + 2a da A 



(3.34) 



where C is the constant appearing in the expression of L2 in (12.490 . In terms of the 
coordinates deflned in (I2.55p . the form Q can be written as: 



2C 



'^10 ^4 



6{p - Pq) dpAVt2A [dy^ A dy^ + dy^ A dy^] 



(3.35) 



Let us now see how we can recover the expression ( I3.35P of Q from a microscopic calcula- 
tion. According to the brane array written at the beginning of this section, the flavor D4- 
branes span a dimension- two surface in the normal bundle A^4. Actually, it is straightforward 
to verify that the family of embeddings ( 12.660 are calibrated by the form K, written in ( 13.130 
and, therefore, they preserve the supersymmetries of the background. The corresponding 
charge distribution obtained by an homogeneous distribution of D4-branes embedded along 
the family (I2.66P and smeared over the [6, (p) two-sphere at a given flxed value pg of the 
transverse coordinate p is given by: 



Kp- PQ)dp A 



sin OdO A. 
47r 



AF 



(3.36) 



where F is the two-form of eqs. (12.680 and (12.720 that results after averaging over the planes 
(I2.66p . Plugging ( I2.72p into ( I3.36P we immediately get (I3.35P if the constant C is related to 
the density nj as: 

C=-2nl,T,^. (3.37) 
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Figure 4: Plots of z{p, a) obtained from the numerical integration of the PDE equation 
f l3.38p for two different values of x = . The left (right) plot corresponds io x = 1/5 
(x = 1/20). In both cases we are taking m pq = 5 and the p < pq region has been obtained 
with the unflavored solution with UV parameters c = — 1 and 2;* = 3. 

By using this result, we can rewrite the PDE (13. lip as: 

2z^z{az + z) = a [z'^ - 2 z z") - a z^'^ 5{p - pq) , (3.38) 
while the discontinuity of ^' at p = pq is given by: 

z'{pQ + e,a) - z\pQ-e,a) = \J z{pQ,(y) . (3.39) 

The PDE fl3.38p with the jump condition fl3.39p can be integrated numerically by using 
the same strategy applied to the case of the D3-branes in section [21 In figure H] we have 
plotted the result of this numerical integration for two values oi nj. From these plots it is 
clear that the effect of the flavors in this solution is qualitatively the same as in the case of 
wrapped D3-branes and, therefore, the remarks made after fl2.75p are also applicable to the 
present case. 

4 A/^ = 2 3d SQCD from wrapped D5-branes 

In this section we will study an alternative way of constructing a gravity dual of three- 
dimensional M = 2 gauge theories, namely by considering the background generated by 
D5-branes wrapping a three-cycle of a Calabi-Yau manifold of complex dimension three. In 
the unflavored case this setup was considered in refs. [HI [15], where an explicit solution 
was found from gauged supergravity. The addition of flavor to this system was addressed 
in ref. [30] by means of a macroscopic approach similar to the one employed here. In the 
present paper we revisit this configuration with our methods and, in particular, we find a 
flavored solution with a good microscopic description in terms of a family of embeddings of 
D5-branes. 

The brane setup we will analyze is represented by the following array: 
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Nf D5 
























In order to formulate a specific ansatz for the metric and RR three-form, let us introduce 
some definitions. Let be SU{2) left-invariant one-forms. We shall consider a rotated 
version of these forms by two angles 6 and (p, with < < vr and < < 27r. Accordingly, 
we define three new one-forms S** {i = 1,2,3) as: 



= sine 
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— sin 


T' 






— cos 6* 
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1 sin<; 
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+ COS( 



w 



w 



S = — COS 9 sin 9 sin 



w 



+ COS( 



w 



Similarly, we define two new one-forms and E'^, obtained from d9 and sin 
the natural frame forms of a two-sphere, as: 



(4.1) 

that are 



= d9 + cos 
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= sin 9 j d(l) 




cos 6' 1 











1 2 

W , W 

h COS (p 

2 ^2 



sm 



(4.2) 



The prototypical example of a non-compact Calabi-Yau threefold with a three-cycle blown 
up is the deformed conifold. It turns out [35] that the metric of the deformed conifold can 
be written in terms of the one-forms and E^ as follows: 



dsl 



^((..)^ + 4(.3)^)+2cosh^(|) ((. 
+ 2sinh^ (I) ({E^f + {E 



l\2 



is 



2\2 



+ 



n2\2 



(4.3) 



where < a < +oo is the radial coordinate, fi is the deformation parameter of the conifold 
and K{a) is the function: 

(sinh(2(T) - 2a)^ 



Kia) 



(4.4) 



23 sinhcr 

We will take the metric (14. 3 p as the starting point to formulate our ansatz for the ten- 
dimensional metric associated to our brane setup. First of all, we add three Minkowski 
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coordinates and a new non-compact coordinate p parametrizing the transverse M 

(—00 < p < +00). As before, we will parametrize the size of the three-cycle by a func- 
tion z = z{p, a) and we will introduce a squashing between the three-cycle and the 5*^ fiber. 
Accordingly, the string frame metric takes the form: 

ds^ = e* (dx^ 2 + ^ dni^ + e-* 1^ [da^ + [{E')^ + {E^] ] + ^p'j , (4.5) 

where $ = $(p, a) is the dilaton and is now: 

= . (4.6) 

In eq. (14. 5 p dVL\ is the metric of the three-sphere, namely: 



dnl 



E (^^)' = \ E (^^)' • (4-7) 



The background is also endowed with an RR three-form F3 which, in complete analogy with 
the results of sections [2] and |3l we will represent in terms of two-forms C2 and A as: 

F^ = dC2 + Q{p~PQ)dk. (4.8) 

We will adopt the following ansatz for C2 and A: 

C2 = giE'AE' + g^S'AS' , g, = g,{p,a) , 

A = LiE^ AE'^ + L2S^ A S\ Li = Li{a) . (4.9) 

Notice that the two-form A introduces a violation of the Bianchi identity of the form: 

dF, = 6{p-pQ)dA, (4.10) 

which corresponds to sources localized at a fixed value pg of the p coordinate. Let us now 
write down the explicit expression for F3. Since we have: 

d{E^ A E^) = -d{S^ A S^) = E^ A AS^ + E^ A AS^ . (4.11) 

Then, the RR three-form F3 given by (14. 8 p and (14.90 takes the form: 

Fs = (^dg, + Q{p-pQ)L,da^ A E' A E^ + {^dg^ + Q{p-pQ)L2da^ A A + 

+ [gi-g2 + e(p-pQ)(Li-L2)) (^e^ as' as"" + e'^as'^as^^ . (4.12) 
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It is natural to consider the following basis of vielbein one-forms: 



e' = e^dx\ (i = 0,l,2), 



eJ = e^^S^-\ (j = 3,4,5) 
m 



z 



e 2 _ E , e = e 2 dp . 



(4.13) 



In order to perform the supersymmetry analysis of the ansatz fl4.5p - fl4.12p in type I IB 
supergravity, let us impose the following set of projections to the Killing spinors: 



r3478 e — e 5 Tsses e — e , Tersg e — —t\ e 



(4.14) 



which leave four supercharges unbroken. In (14.141) t\ is the first Pauli matrix. The corre- 
sponding BPS equations are: 
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3m 



e^* ^ ( fi'i + 0(P - Pq) i^i ) + ^ (^71 - + e(p - pq)(Li - L2)) 



i = 2m2e-2*- 
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z' = 2m^^{g2 -gi + Q{p - Pq)[L2 - U)) 



92 = --{gi - 92) - 0(p - pq] 
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L2 — Li 
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-2<S> 



a 



92 = -r 

'z 



The consistency of this system requires that the Lj's satisfy the condition: 

Li — Li 



a 



(4.15) 



(4.16) 



Moreover, from the BPS system f l4.15p we get the following PDE for the function z{p, a): 



z \ z z + -z" 



^12 II . (L2 — Li) 0/2 ^, , 

-z'^- z z" + ^ '- z""^ 6{p- Pq) 

I cr 



(4.17) 



One can check that, if eqs. f l4.15p and f l4.16p hold, then the equations of motion of the gravity 
plus smeared branes system are fulfilled. Notice the difference between the consistency 
conditions f l4.16p and (12.471) . which is due to the fact that in this D5-brane case the flavor 
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branes are co dimension one in the normal bundle N^, whereas in the setups of sections [2] 
and [3] the flavor branes span a codimension two submanifold of A^4. 
The generalized calibration form for this case is a six- form of the type: 

JC = ^/C,,...„,e'^i-"« . (4.18) 
The different components of IC are given by: 

^ai-as = ...ag Ti e , (4.19) 

where e is a Killing spinor of the background, normalized as e~ 2" e' e = 1, and the minus sign 
in the definition fl4.19p has been introduced for convenience. By using the SUSY projections 
satisfied by our solutions (eq. f l4.14p ). we get the actual components of K, in the frame basis 
fl4.13p . namely: 

JC = -e°12 ^ ^^345 ^ g367 ^ g468 _ ^578) _ ^4 ^0) 

The BPS system 04.151) can be now summarized in the following equations satisfied by the 
calibration form /C: 

*F-i = -d(e-*/C) , d(*/C) = . (4.21) 



4.1 The unflavored solution 

The background of refs. [HI [15], obtained in seven- dimensional gauged supergravity and 
afterwards uplifted to ten dimensions, provides a solution of the BPS system fl4.15p and of 
the PDE equation fl4.17p in the unflavored case Li = = 0. Let us present this solution, 
entirely written in our variables. The function z{p, a) is given as the solution of the implicit 
equation: 



x3/2 (J, 



+ 



^i^3/4(f)) v^(/-l/4(f)+ciri/4(f))' 



4m^ 



(4.22) 



where /,, and K,, are modified Bessel functions. One can, indeed, check that the function 



z{p, a) defined by (141221) solves the PDE equation f l4Tfl) for Li = La 
the function x(z) as: 



1-1/4(1) + cK 



'3/4I2, 



0. Let us next define 



(4.23) 



Then, the dilaton $(p, a) is given by: 



-6x n2 



-1 q2x gr2 



(4.24) 
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In order to write down the functions gi and g2 that parametrize the three form F3, let us 
now define the angle ip [0 < ip < tt) as: 

cotV = z5e-2^^ . (4.25) 
a 

In terms of ^p the functions gi and (72 are simply given by: 

^i = ^- e ^2 = ^- 4.26 



This unflavored solution was analyzed in detail in ref. [15], where it was argued that, for 
c > 0, it has a good IR singularity corresponding to a linear distribution of fivebranes. By 
means of a probe computation, it was proposed in [I5] that the solution describes a slice of 
the Coulomb branch of A/" = 2 , = 3 pure YM theory. Actually, one can also prove in this 
case that the SUSY locus of color D5-brane probe occurs at cr = and that the Yang-Mills 
coupling is related to z{p, a = 0) by: 

9ym (An)'^ gs a' 

Notice that now z does not approach a constant value in the UV. Indeed, from f l4.22p one 
gets that it grows logarithmically as z{p,0) ~ 21og(mp) for large p. Accordingly, gYM in 
f l4.27p decreases as a power of a logarithm of p in the UV, contrary to the result expected 
for a three-dimensional gauge theory. Notice also that the dilaton blows up in the UV and, 
therefore, a string theory completion is needed in this region. 

4.2 Flavored solutions 

The modification of the Bianchi identity of F3 induced by a set of flavor D5-branes distributed 
with a charge density four-form Q is: 

dFs = 2Kl^T^Vt . (4.28) 

By comparing f l4.28p with the value of ciFs for our ansatz (eq. f l4.10p ). one immediately gets 
the smearing form VL in terms of A, namely: 

n = ^(P - Pq) dp A _ ^^ 29) 

As in the D3- and D4-brane cases studied above, we will concentrate on studying the config- 
uration in which the function Li that parametrizes A vanishes. In this case the compatibility 
condition fl4.16p can be easily solved: 

Li = 0, L2 = Aa, (4.30) 
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with A being a constant. The smearing form for this configuration can be straightforwardly 
obtained by using (14. lip to compute dA. One gets: 

n = ^4^PfI^dp A{daAS'AS^-a {S' A A E' + A A E')) . (4.31) 

If we now substitute in (14.311) the expression for the one-forms S*^'^'^ and E^'"^ written in eqs. 
(14. ip and (14. 2p . we see that they can be rearranged so that Q can be recast as: 

n = dpA [d (cr sin 9 sin (p) A w"^ A + d {o sin 9 cos (j)) Aw^ A 

+ d{acos9) Aw^ Aw"^]. (4.32) 

Let us introduce now the following set of cartesian coordinates for N^: 

= (J sin 6^ sin , y'^ = cr sin 6^ cos , y^ = acos9 . (4.33) 

In terms of these cartesian coordinates the smearing form has the following appealing and 
compact form: 

Q = -^^iP^lRl e^^k^p^dy^ ^w'' . (4.34) 
16 -is 

In the next subsection we will identify a set of supersymmetric embeddings which give rise 
to a charge density distribution given by (14.340 . 

4.3 Microscopic analysis 

First of all, let us determine a continuous set of supersymmetric D5-brane embeddings dis- 
tributed as the flavor branes of our setup. In these embeddings the D5-brane is located at a 
fixed p = pq, is extended along a two-dimensional plane in N^, and, simultaneously, wraps an 
5*^ inside the three-cycle S^. The particular 5*^ C that the brane wraps depends on the 
plane that it spans in N^. In order to describe this relation, we first define the three- vector 
n as: 

n = (sin a sin /3, sin a cos/3, cos a) , (4.35) 

where < a < vr and < /3 < 27r. Let us next consider the family of embeddings in N3 
defined by the following function: 

f(y)=n-y-y, = 0, (4.36) 

where y = {y^,y'^,y^) are the cartesian coordinates defined in (14.330 . This is the equation 
of a plane having n as its normal direction and with y^, being its minimal distance to the 
y = origin of N^. Let us next define the following two tangent vectors to the plane: 

ti = (cos /3 , — sin /3 , 0) , 

t2 = (cos a sin /3, cos a cos /3, — sin a) . (4.37) 
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Notice that {ti , t2,'n) is an orthonormal basis in R^. Let us now arrange the three SU{2) 
one-forms in a vector, namely: w = {w^ , w'^ , w"^) . If we define the one-forms: 

wt, = U ■ w , Wn = n-w , (4.38) 

then, the embeddings in the are defined by the two differential conditions: 

wt, = , wt, = . (4.39) 

Notice that the system fl4.39p is integrable due to the property: 

dwt^ = -wt^ A Wn , dwt^ = -Wn A Wti , (4.40) 

which shows that the differential equations derived from fl4.39p are on involution and, ac- 
cording to Frobenius' theorem, they are integrable. Moreover, in order to show that the 
corresponding embedding preserves the supersymmetry of the background, one has to verify 
the fulfillment of the generalized calibration condition: 

/C = Vo\{Mg) , (4.41) 

where /C is the pullback of the six-form f l4.20p . To verify fl4.4ip we will come back for a while 
to our original (cr, 6, (p) coordinates of A^s. Notice, first of all, that we can use the differential 
equations fl4.39p to express the pullback of two of the w*'s in terms of the third one. For 
example, if a 7^ 7r/2 we can write: 

vj-^ = sin/3 tana w^ , = cos/3 tana w^ . (4.42) 

Using this result we can compute the pullback of the different one-forms of the basis fl4.13p 
and, thus, one can compute the pullback of the calibration form /C. The result can be written 
as: 

^ ^ a e'f' sec a sm 9 ^ ^ ^^2 ^'^ dO A d(j) . (4.43) 

2m^ \/2z cos a cos 9 + cos(/3 — 0) sin a sin 9 

One can also compute the induced metric and check that, indeed, eq. f l4.4ip is satisfied and, 
therefore, the embeddings described above are supersymmetric, as claimed. 

Let us now verify that the charge distribution generated by this family of D5-brane config- 
urations is just given by a six-form Q such as the one written in (14. 34 p . Following the same 
methodology that we applied for the D3- and D4-brane cases, let us write the smearing form 
Q as: 

n = 5ip- pq) dp at, (4.44) 
where F is the following three-form: 

F = / dp6{f) [dfATs^], (4.45) 



34 



with the integration measure djj, given by: 

sin a da df3 

dfi = rifdy^ , (4.46) 

and with F^s being the charge density distribution obtained after averaging over the possible 
embbedings on the for a given plane in N3. In (14.46P n/ is a density of flavor branes. The 
integral over in (14.451) can be immediately performed by using S{f) = 5{n ■ y — y^). One 
gets: 



Att 



j sin a da d(3 [df ATgs] . (4.47) 



In order to determine Ygs let us first consider the particular plane with a = f3 = 0. In 
this case the differential equations for the embedding in the are just = w'^ = 0. After 
looking at the parametrization ( ]2.54p of w-^ and w'^ in terms of the three angles {6, 0, one 
immediately realizes that the equation for this embedding can be integrated as: 

/.^ = - 0, = , fg = 9 - 9, = 0, (4.48) 

where 0* and are constant angles. These embeddings depend on two continuous parame- 
ters (j)^, and 6^, which span a two-sphere. Therefore the corresponding two-form F^s is given 

^^^(^*d<^*d(^*SiU)6ife) [dUAdfe] , (4.49) 



S3 



Air 



where we are integrating over 0* and 6^ with the invariant measure of the two-sphere. These 
integrals are immediately performed with the help of the Dirac ^-functions: 

sin B ^ 1 

F53 = —— d(j)Ad9 = —w^Aw\ (4.50) 
Ati An 

From this result it is clear that the generalization to arbitrary values of a and /3 is: 

F53 = -^u'ti Awt2 . (4.51) 
An 

Taking into account that n = ti x t2 or, in components, = e^^^t\t2i we can write: 

Wt^ Awt^ = ^e'^'^n'w^ A . (4.52) 
Thus, we can rewrite (I4.5ip as: 

F53 = ^e'^^n'w^ A , (4.53) 
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and, plugging this result in f l4.47p . we find that T can be represented as: 



r = ^ 

327r2 

Moreover, using the fact that: 



27r 



sin a da d/3 



(4.54) 



sin a da 



2n 



dl3 n' = — S'^ 



we get: 



e^^^ dy' A A 



247r 



(4.55) 



(4.56) 



Therefore, by combining this result with f l4.44p we conclude that the smearing form Q for 
this family of embeddings can be written as: 



n 



f Ajk 



247r 



5{p - pq) dp A dy' A A 



(4.57) 



By comparing fl4.57p and the macroscopic expression of written in fl4.34p we get that both 
expressions coincide if the constant A is identified with: 



A 



"2^10 ^5 



37r 



(4.5^ 



Taking into account that 2k^q T5 = Ati"^ a', we can rewrite the previous relation as: 



A = — — gsa'uf . 



(4.59) 



4.4 Numerical results 

Using the expression fl4.59p of the constant A, one can straightforwardly demonstrate that 
the PDE fl4.17p for the flavored system takes the form: 

. ^ + i i^) = i - . - ^ z'/^ 6{p - pq) . (4.60) 

The discontinuity of z' at p = pq in this case is given by: 



z\pq + e, (t) - z\pQ - e, d) = \J z{pQ^ ^) • (4-61) 

In figure [5] we show the numerical result of the integration of the PDE fl4.60p for z{p, a) 
which reduces to the unflavored solution for p < pq and such that its first partial derivative 
with respect to p has the discontinuity given in (14.610 . 
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Figure 5: Plots of z{p, a) obtained from the numerical integration of the PDE f l4.60p for two 
values of a; = The left curve corresponds to x = 2/3, while the curve on the right was 

obtained for x = 1/5. In these two cases m pq = 3 and we are considering the case with 
c = 1 for p < Pq, where c is the constant appearing in the unflavored solution fl4.22p . 

5 Summary and conclusions 

In this paper we have studied supergravity backgrounds which are obtained by wrapping 
D-branes on cycles of a Calabi-Yau cone. In their unflavored version these backgrounds were 
obtained some time ago by uplifting to ten dimensions the corresponding solutions in lower 
dimensional gauged supergravity. These uplifted solutions depend on several angles and are 
quite involved. Here we have rewritten them in terms of the two radial coordinates (p, a). In 
spite of the fact that one has to deal with functions depending on these two variables, this 
form of presenting these backgrounds is a simplification that clarifies their interpretation and 
allows to easily find their generalization when flavor sources are added. In the cases of the 
D3- and D4-branes, the solutions of the BPS systems in gauged supergravity of ref. [5] can 
be obtained by solving the master equations fl2.18p and f l3.15p and using these results in the 
implicit equations (I2.19P and (I3.16P respectively. In the case of the D5-brane, the analytic 
solution of the unflavored background of refs. p31 [15] is obtained, in terms of our variables, 
by solving the implicit equation f l4.22p . 

The most important results of this paper are the construction of backgrounds which in- 
clude the backreaction due to a large number of flavor branes. The latter act as extended 
dynamical sources for the different supergravity fields and, in particular, they modify the 
Bianchi identity of the RR fields. In order to deal with this problem, we have first adopted 
an effective macroscopic approach in which the precise knowledge of the flavor brane embed- 
dings is not needed. 

Using the preservation of supersymmetry and the compatibility with the metric ansatz as 
guiding principles, we were able to parametrize the violation of the Bianchi identity induced 
by the flavor branes by means of two functions Li and L2 and we argued, as in ref. [IT], that 
the most physically sensible configuration corresponds to taking Li = 0. We have verified 
this fact by reproducing the results found in the macroscopic approach by means of a careful 
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analysis of the supersymmetric configurations of the flavor branes. We first found a family of 
embeddings which satisfy a generalized calibration condition and such that, in an appropriate 
system of coordinates, they can be regarded as planes in the directions orthogonal to the cycle 
wrapped by the color branes. We then computed the RR charge density fl resulting after 
distributing homogeneously the flavor branes throughout the family of calibrated planes and 
we checked that the result coincides with the one found in the macroscopic approach. In turn, 
this result allowed us to relate the constants A and C of the macroscopic approach to the 
density ri/ of flavor branes. Once this identification is performed, the flavored backgrounds 
are obtained by integrating numerically the PDE equation for z{p, a) with the condition 
that the solution reduces to the unflavored one for p < pg and such that the first derivative 
z'{p,a) jumps at the location p = pq of the flavor brane as dictated by the corresponding 
BPS equation. 

Let us now comment on some points not addressed in our work that could be interesting 
to analyze in the future. First of all, we could study the fluctuations of a flavor probe brane 
embedded along one of the supersymmetric planes. By studying the normalizable modes 
one would find the meson mass spectrum. The results should be close to the ones analyzed 
in refs. [231 [IHl EI] for other backgrounds similar to the ones analyzed here. In particular, 
by comparing the spectra for the quenched and unquenched backgrounds one could study 
the effects of the backreaction on the meson masses. When the probe brane is located at 
p = Pq < Pq the spectrum for the backreacted solution is the same as the one without 
backreaction since both metrics coincide for p < pq. On the contrary, for pq > pq this is 
not the case and, from the change in the mass spectra, one can extract, in the holographic 
setup, the corrections produced by charge screening effects due to quark loops. As in other 
similar cases analyzed before ^31 [161 EU [29], one expects these screening effects to be small. 

In this paper we have only studied configurations of flavor branes extended along cal- 
ibrated planes. There exist, however, more general embeddings that preserve the same 
supersymmetry (see, for example, section [2.4.1l ). Generically, these other embeddings repre- 
sent configurations in which the flavor branes end on the color branes in such a way that both 
types of branes are recombined. On the field theory side these solutions should correspond 
to Higgs branches of the corresponding supersymmetric theory. 

One aspect that is missing for the kind of gravity duals analyzed here is the precise dic- 
tionary between them and the corresponding field theory. If we had this dictionary we could 
extract more precise information from our solutions. We could, for example, understand the 
relation between the field theory duals to the backgrounds of sections [3] and [H which should 
correspond to two three-dimensional gauge theories with different UV completions. Nev- 
ertheless, we regard our work as part of the efforts devoted to approach the gauge/gravity 
duality to particle physics phenomenology. In the absence of a true gravity dual of QCD, the 
best that one can do is trying to uncover universal features of the duality by systematically 
analyzing different backgrounds and their generalizations. We hope that some of the results 
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presented here could be useful in this endeavor. 
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